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ABSTRACT 


This thesis was undertaken to determine the effects of axial 
and radial clearances on the performance of a single stage turbine 
with blunt leading edges and non-twisted blades. A series of tests 
was conducted on the so-called Mod II Turbine using the Transonic 
Turbine Test Rig of the Turbo-Propulsion Laboratory, Department of 
Aeronautics, of the Naval Postgraduate School. The results of these 
tests are presented together with a comparison of the experimental 
results and results predicted by a three-dimensional turbine perform- 
ance calculating method. In addition, measured flow conditions up- 
stream of the stator, between the stator and the rotor, and at the 


rotor discharge are presented and compared with predicted values. 
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l. Introduction 

The Mod II Turbine is a test model of the fuel pump drive turbine 
designed by Professor M. H. Vavra of the Department of Aeronautics, 
Naval Postgraduate School, for an advanced liquid rocket engine. It 
was designed for good off-design point performance, ruggedness, and 
simplicity. Good off-design point performance was needed because the 
rocket was to produce variable thrust over a large percent variation 
of full thrust. Ruggedness was a design criteria because the rocket 
must endure stop-start operations. Simplicity of design provides ease 
of manufacture and low unit cost. In addition, the Mod II Turbine was 
designed to allow for adequate space in the blades for internal cool- 
ing passages. The foregoing considerations lead to a turbine design 
with thick non-twisted blades which have blunt leading edges. 

The amount of test data which is available on turbine designs of 
this type is very limited. Turbine designers have a particular need 
to know what effects are induced in turbine performance by varied 
clearances between the stator and rotor and between the rotor tips 
and the shroud. This project hopefully helps fill this need, 

The Transonic Turbine Test Rig located at the Aeronautics 
Propulsion Laboratory of the Naval Postgraduate School provides a 
unique test bed for determining the effects of axial and radial clear- 
ances on the performance of turbines at varied pressure ratios and 
speeds. This report covers the testing of the Mod II Turbine in the 
Test Rig at two radial clearances, five axial clearances, four pres- 
sure ratios, and varied rotational speeds between 10,000 and 19,000 
RPM. 

The author is very appreciative of the many hours of cheerful 


assistance during the equipment set-up and data acquisition given by 
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Mr. Jim Hammer of the Department of Aeronautics. Also, the willing 
assistance of Lt. P. M. Commons, U.S.N. and the helpful advice given 
by Lt. R. H. Harrison, U.S.N. are greatly appreciated. Particular 
thanks are given to Professor Vavra for his patient instruction and 


guidance. 


2. Turbine Description 


Figures 1, 2, and 3, and Figures 4, 5, and 6 are photographs of 
the Mod II Turbine stator and rotor, respectively. The turbine is a 
single stage unit of the reaction type. Pertinent turbine dimensions 
are listed in Table I together with the design point performance 
parameters. Unless otherwise specified, the values in Table I refer 
to the mean blade radius. Most of the dimensions given in Table I 
are self-explanatory; however, the following quantities need to be 
defined more clearly. The throat diameter, a, shown in Fig. 7, is the 
minimum distance between the blades at the mean radius. The throat 


area, Avy is defined as 


= (1) 
Aun = Zah 


Definitions of the absolute and relative exit angles, & and g , are 
shown in Fig. 8. Figure 8 depicts the so-called velocity triangles 
for the Mod II Turbine at the design conditions. All angles are 
measured from the axis. Positive angles are measured to velocity 
vectors with peripheral components in the direction of the rotor 
rotational speed vector, ill Definitions of the loss coefficients 

and other performance parameters listed in Table I are presented in 
Section 4, During the designing of a turbine, loss coefficients must 


be assumed in order to determine proper areas for the flow channels. 
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Large losses can result from improper area design, Turbine designers 
usually use experimentally determined loss coefficients. The loss 
coefficients used in designing the Mod II Turbine and listed in Table I 


were derived from experimental data presented by rein 


3. Test Installation 

The turbine was tested on the Transonic Turbine Test Rig, here- 
after referred to as the TTR. The TIR is designed and instrumented to 
determine the performance characteristics of a turbine stage at varied 
axial and radial clearances, and to measure the flow properties before 
and after each blade row. It has the unique capability of determining 
mean flow conditions in a turbine without introducing flow disturbing 
probes into the flow stream. Detailed descriptions of the TTR instal- 
lation, instrumentation and data reduction techniques are presented by 
egemonsce This discussion will be limited to the salient features of 
the TTR from which performance data on the Mod II Turbine were obtained. 

Compressed air is used as the working fluid in the TTR. It is 
supplied by a twelve stage axial compressor located in a test cell 
adjacent to the TTR. The supply air enters the TIR test cell at the 
inlet valve for Tank 1 as shown in Fig. 9. Tank 1 is a settling tank 
and plenum chamber for the exhauster and for Tank 2. The exhauster is 
used to lower the pressure in the test hood below atmosphere to obtain 
large pressure ratios across the turbine being tested. 
SKalener Armin, Experimentelle Nachprtfung eines Berechnungsver- 


fahrens fur axiale Stromungsmaschinen am Beispiel einer Turbinenstufe. 
Forschg. Ing-Wes. 31 (1965) Nr. 5. 


commons . P. M., Instrumentation of the Transonic Turbine Test Rig to 
Determine the Performance of Turbine Inlet Guide Vanes by the Application 
of the Momentum and Moment of Momentum Equations (NPGS Thesis, Sept. 1967). 
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Air for the turbine leaves Tank 1 through the flow rate nozzle 
and goes into Tank 2 which is the settling tank for the turbine plenum. 
The flow rate nozzle is instrumented with total pressure taps and a 
total temperature probe ahead of the nozzle and static pressure taps 
downstream of the nozzle. The air path from Tank 2, through the 
turbine inlet valve, and into the test hood to the turbine plenum is 
shown in Fig. 9. 

Figure 10 is a cross-section of the arrangement inside the test 
hood. The turbine plenum is supported independently from the stator 
assembly and surrounds the upstream end of the stator assembly. The 
plenum is instrumented with total temperature and total pressure probes. 
A small amount of air leaks from the plenum through the labyrinths that 
seal the plenum presSure from the pressure inside the hood. Air flows 
radially from the plenum into the stator assembly, thence through the 
turbulence reduction screens to the stator. Six fixed total pressure 
probes and two moveable total pressure and temperature Kiel probes are 
located between the screens and the stator. The total pressure probes 
are fixed at a radial distance which approximately divides the flow 
entering the stator into equal mass flow rate increments. The fixed 
probes and the conical turbulence reduction screens are shown in 
Rage, ll. 

The closure plate, shown in Fig. 12, is attached to the stator 
assembly by a cylindrical member and a spoked wheel type flexure 
device, which is instrumented with strain gages to determine the 
torque and the axial force transmitted from the closure plate to 
the stator assembly. Also, the cavity between the closure plate 


and the stator assembly is instrumented with a static pressure port. 
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Static pressure ports are arranged along the shroud at %-inch 
intervals beginning at the exit of the stator. The arrangement of 
these ports on the downstream end of the shroud is shown in Fig. 12. 

The radial clearance between the rotor tips and the shroud is 
determined by the inside diameter of the shroud. Different shrouds 
are available for installation on the test rig to provide varied 
radial clearances. 

The entire stator assembly is supported by flexures which allow 
the assembly to move in the axial direction and around the axis. 
Reluctance type force gages are located between the moveable stator 
assembly and the fixed structure. The force gages measure the axial 
force, and the moment about the axis, which are transmitted from the 
moveable stator assembly to the fixed structure. 

The rotor is cantilevered from the rotor bearing support stand. 
Figure 14 shows the rotor mounted in the bearing stand. Two sets of 
matched precision ball bearings support the rotor shaft. The bearings 
are lubricated by an oil mist system. 

The axial distance between the stator and rotor is varied by 
sliding the rotor bearing assembly in the support stand. The bolt 
at the top of the rotor bearing support stand locks the sliding 
assembly in the desired axial position. 

A quill shaft connects the rotor shaft to the dynamometer shaft. 
A six-spoked flux cutter is fitted to the dynamometer end of the quill 
shaft. The flux cutter passes through the field of a magnetic pickup 
from which the RPM of the rotor is obtained. 

The dynamometer is an air brake device which is cantilevered 
from the dynamometer bearing support stand. The dynamometer bearings 


are similar to the rotor bearings. A twenty inch long arm is attached 
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to the shaft between the dynamometer bearing stand and the dynamometer. 
The arm acts on a reluctance type force gage from which the dynamometer 
torque is determined. Table II lists the quantities recorded during 

a performance test run. All pressures are measured on a mercury 
manometer, except the pressure differential across the flow nozzle 
which is read on a water filled, U-tube. All the temperatures are 


obtained from Iron-Constantan thermocouples referenced to an ice bath. 


4, Analysis and Data Reduction. 


General 

The turbine performance data were analyzed with a one-dimensional 
mean streamline approach. Steady axisymmetric flow conditions were 
assumed to exist at the entrance and discharge of each blade row. 
Adiabatic flow conditions were assumed to exist through the entire 
stage. The mean flow conditions were assumed to exist at the root 
mean square of the blade radii, hereafter referred to as the mean 
radius, R The mean radius is found as 


2 2 | 
Rin = (A SF — inches) (2) 


where Ry = radius of blade tip (inches) 


RL = radius of blade hub (inches). 
For simplicity in the following development, the mean radius will be 
considered constant through the stage. In the actual calculations, 
the variation in mean radius between the stator blades and the rotor 
blades was taken into account. The major portion of the data reduction 
calculations Was performed by the IBM 360 digital computer located at 


the Naval Postgraduate School, 
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Flow rate 
The flow rate through the turbine is the flow rate through the 
flow measuring nozzle less the flow rate through the labyrinth seals. 
Thus 
i = Wy - W, (3) 
where 


W = Flow rate through the test turbine (1b_/sec) 


Se 
ul 


Flow rate through the flow measuring nozzle (1b_/sec) 


L Flow rate through the labyrinth seals (1b_/sec) 


The flow rates through the flow nozzle and through the labyrinth seals 
are determined by measuring the applicable quantities listed in 
Table II. The calibration methods and data reduction formulae are 
discussed by Comons 
Stage entrance properties 

The total pressure at the stator entrance, Po? is taken as the 
average pressure indicated by the six fixed pressure oPenee.. The 
validity of this assumption is explored in Sec. 7. The total temper- 


ature, T. , is obtained from the Kiel probes. 


jee. 
Stator discharge properties 

The fluid properties at the stator discharge and the stator 
performance parameters are obtained by determining the velocity triangle 
at the stator exit. The peripheral component, Muy? of the absolute 
velocity is obtained by application of the moment of momentum law to 
the fluid contained in the stator assembly. From the derivation 
presented in Appendix I, 


3 
ibid, Sec, 4 
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Vu, = (Ma +M5)9 AWRm (4) 


where 


V = Peripheral component of absolute velocity at the stator 
u 
i 
exit. (ft/sec) 


M. = Moment applied to the stator assembly by the moment capsule 
shown on Fig. 12. (f£t-1b,) 


Mor, = Moment applied to the stator assembly by the closure plate 


(f£t-1b_) 
The axial velocity component, ve , 1s found by application of 
if 


the conservation of momentum law to the fluid in the stator assembly. 


From Appendix I, 


ke 
a ( ] (5) 
- = r 
V2 ale th +2h- 20 Pad 
Ih, 
where 
ve = Axial velocity component at the stator discharge (ft/sec) 

1 


By = Force applied to the stator assembly by the axial force gage 
(See Fig. 13) (1b ,) 
For. = Axial force applied to the stator assembly by the closure 
plate (1b ,) 
F, = Net pressure force acting on stator assembly. (1b -) 
(See Appendix I, p.112). . 
p = Static pressure at any radius at the stator exit. (1b ,/in*) 


Refer to Appendix I, p.112 , for the solution to the pressure integral 


given in Eq. (5). 


Zu 


Since the continuity equation must be satisfied by the fluid flow 


through the stator, the axial velocity component, vee , may be determined 
1 
as, 





: sep At PY _2we, Wi Ye 
Ma oleh ee Cl+ (eta - WTP 


Refer to Appendix I, p.112 , for the derivation of Eq. (6). 

In the data reduction computer program, contained in Appendix II, 
both Eq. (5) and Eq. (6) are satisfied by an iteration procedure which 
determines a compatible static pressure distribution in the radial 
direction at the stator exit. The stator exit velocity is then found 


as 


l 
V, = [ Vas emi? 7) 


Referring to Fig. 16, the stator absolute discharge angle is 
: Vi 
= Tan! (i) (8) 


The relative axial velocity component is 


Va, (9) 


and the relative peripheral component is 


e - (10) 
where 
@ = rotor rotational speed (rad/sec) 
The relative velocity is now found as 
= [ W, +w?]? (11) 
a Vy 


and the relative stator discharge angle is 


B, = Tane! ( ) (12) 


Aj 
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The absolute and relative Mach numbers are found as 


V4 13 

M, = ry Glsy, 

- W (14) 
me Qy 


where a, is the local speed of sound at the stator exit, obtained as 


il 
My, 
QA, = C98R. Ty, ] Ar 


where 


¥ = specific heat ratio for air (1.401). 


The effectiveness of the stator is measured by the stator loss 


coefficient which, referring to Fig. 15, is defined as 


ie Ti BR cs This (16) 
Ss Tee Treg 





where 
= static temperature at the stator discharge if the expansion 


Me 
is 
iE 1 EC. 
process from (To ee to p, were isentropi 


It can be seen on Fig. 15 that an equivalent definition for the 


Stator loss coefficient is 


2 2 
C _ Mun (17) 
a 2 
Ss Ving 
where 
Vi velocity at stator exit if expansion through the stator were 
th 


isentropic, 
By using the energy equation and the equation of state for an isentropic 


expansion 


f=! 
Vy, =2gdep ko [1 (Mp d* J (18) 


Zo 


Combining Eqs. (16) and (17) 
- 1 Me 

3 I Gate ICR DFT 09 
¢, is an average value of the stator effectiveness. In the actual 
flow passage through the stator the major part of the losses is due 
to the boundary layers which form along the surfaces of the blades. 
In fact, expansions along streamlines away from the bounding surfaces 
should be nearly isentropic. For this reason, it is quite difficult 
to obtain a meaningful average stator loss coefficient by measuring 
the flow properties with probes introduced into the flow stream. 
Obviously, it would require a great number of points to be measured 
very near the surfaces. Herein lies the significant advantage of the 
TTR with the force capsule arrangement which gives mean flow conditions 
by measuring external forces on the stator assembly. The stator 


efficiency is given by 


Von LS $ (20) 


The flow function d is given by vewee as 


bef 
b ; re Ath Ro = y" soi 


and the isentropic flow function as 


y 
6 =f 2h (KO -(IP we 


The stator flow restriction factor or blockage factor is defined as 


B= ¢ (23) 


VetE ey M. H., Problems of Fluid Mechanics in Radial Turbomachines 
Parts I & II. Von Karman Institute Course Note 55a. Rhode-Saint- 
Genese, Belgium: Von Karman Institute for Fluid Dynamics, March 1965, 
equation C(7). | 
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Rotor discharge properties 
By considering the flow through the rotor relative to the moving 
rotor blades, the fluid properties at the rotor exit can be determined 
in a manner quite analogous to the method applied to the flow through 
the stator. 
The peripheral component of the relative velocity at the rotor 
exit is found from the moment of momentum equation, From Appendix I, 


peeee2, © 


Wy, = Wy - Mo ey (24) 
m™ 


where 
My = moment measured by the dynamometer force capsule. From 


Appendix I, pll2 , the axial component of the relative velocity is 


WRe Te 
VV aes SiS (25) 
Qo P Ay 


where the static temperature at the rotor discharge is found by 


combining the energy and continuity equations as 


Tx Weep Ar Py i Depa ft 4 _2W? Re 


2° W? Re Wis 24h _Wi_y] “4 (26) 


ap At PE aa dep 24 J<p re 


The various velocity components and discharge angles for the rotor 


exit can now be computed. From Fig. 16, 


Y, 
Wome (Wa, +W,, : 27) 
Va. = Wa. (28) 


INO 
Wi 


¥, 
[Vere Viel (30) 


V. = 
Wy 
6, = Tage a (31) 
oA 
o, = TAN (H) (32) 
2 


Analogous to the stator loss coefficient, the rotor loss coefficient 
is defined as 


C 4- Wa 
R We, (33) 


W, = relative velocity at the rotor exit if the expansion 
th 
through the rotor were isentropic, (ft/sec). 
This expansion is shown in Fig. 15 from (P37) to (PT, » 
{ is 


W, is calculated as 
th 


: y, 
Za ie in ig 
We, = {mit e2gJeqp[1- (39) IF ae 
The efficiency of the rotor is defined by 


Ne =p oe (35) 


Stage Performance Parameters 
The overall efficiency of the turbine stage is that percent of 
the isentropic temperature drop across the stage which is used in 
developing work output. The isentropic temperature drop across the 


stage, at , is shown in Fig. 15. It is computed as 
8 
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a 
ATa = Tell (%)" 1 (36) 


AT. . is often expressed in terms of a theoretical velocity, 
C ,as 
2 
Co 
f= == (37) 
At 29J &p 
The temperature drop which represents the work output is shown on 


Pre, 15 as AT . It is computed as 


AT, = — (38) 


Or, bY Eulers turbine equation, as 


a ie U (V1 -Vun) GJ ep 2 


The efficiency is then computed as 





Aly .. £U Vu~Wv2) Mp W/W &p 

nN = ~ = vy (40) 
AT : he L1-C%/) **] 

to 


This efficiency is known as the total to static efficiency. 
The isentropic expansion across the stage is considered to start at 
the total pressure ahead of the stage and extend to the static 
pressure after the stage. The kinetic energy of the fluid leaving 
the stage is considered as a loss. 

The percent of AT... associated with a theoretical isentropic 
expansion across the rotor from Py to P» is known as the theoretical 


degree of reaction. In terms of velocities it is defined as 


2 


4, 
Oh ( 1) 


— 
% 
{1 


27 


vw 
r is an important stage parameter. It is an indication of the 


amount of acceleration experienced by the flow in each row of blades 
‘ 
in the turbine stage. Using Eqs. (18), (36), and (37), i is expressed 


as 


vy 
re Gap = 1 (42) 


(Ryp ye - 1 
The values of na at the hub and tip are computed by substituting Ph 
and Py for Py> in Eq. (42). 

The peripheral speed of a turbine rotor, U, is usually a fixed 
value which is determined by the allowable stress level in the rotor 
blades. For this reason, U is commonly used to define dimensionless 
stage performance parameters, 

The isentropic head coefficient, kg? relates the isentropic 


energy change across the stage to the peripheral speed by the relation 


Co \% 
Rip= Gee (43) 
The work coefficient k relates the actual work accomplished in 


a stage to the peripheral speed as 


_ 2(W,-W2) ~  =«_Alw 
Rag ee (44) 


U VAgT cp 
Keg is used by designers to estimate the number of stages 
necessary to handle a given isentropic energy change at any given 
peripheral speed, U. Similarly, kK is used to estimate the flow rate 
through a stage necessary to produce a specified amount of power at any 


given speed U, 
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In order to compare the results of the performance tests from run 
to run, and with predicted results, the NASA reference system is 


employed. The reference parameters are defined as 


Giimas”™* eK 401)(53-35)C 518-4) (45) 
fOr alr 
- dee 
Q@ = "518.4 (46) 
and 
Re 
6 = “Yat (47) 


The performance values which are measured on the TTR and then 


reduced to referred values are flow rate, 


Waa W 1g (1b_/sec) (48) 


dynamometer moment, 


Moses Mo Ca (49) 
horsepower, 
HPo-. = oe (hp) (50) 


and rotational speed 


Neer = (rpm) (51) 


“0 

5. Description of Performance Tests 
The performance tests were conducted from January through August 

1967. The turbine was operated about 150 hours while performance data 


were being taken. Table III lists the test parameters used in each run, 


The first series of test runs was performed at a radial clearance 
between the rotor tips and the shroud, AP , of 0.033 inch. Prior 

to each run, the axial clearance between the stator and rotor, Ax, 
was set at the desired position. Five axial positions, Ax, of 
0,200, 0.0410, 0.620, 1.000 and 1.500 inches, were used in the first 
set of tests. The second series of tests was performed at a radial 
clearance, AF , of 0.015 inch. Since the changes in performance 
parameters with a change in axial clearance were found to be small 
at the original radial clearance, it was decided to conduct the final 
set of tests at only three axial clearances} namely, Ax of 0.200, 
0,405. andar, 000. inch, 

At each combination of axial clearance, Ax, and radial clearance, 
Ar, the turbine was tested at four pressure ratios across the stage, 
(Pip) namely, at 1.3, 1.4, 1.5, and 1.6. An additional pressure 
ratio, (PB 4 /P) of 1.45, was examined during run number 63. 

The test data items shown in Table II were recorded at a number 
of rotational speeds, N, between 10,000 RPM and 19,000 RPM, at each 
combination of radial clearance, axial clearance and pressure ratio. 
The pressure ratios listed in Table III are approximate due to the 
difficulty experienced in holding exact pressure ratios while varying 
the RPM. This was a particular problem during the initial tests 
conducted with the exhauster operating. As experience was gained and 
better techniques were developed, it became possible to hold the pressure 
ratios to within one-half a percent of the desired pressure ratio. 

The torque absorption capacity of the dynamometer limited the 


range of head coefficients, k over which the tests could be performed. 


is’ 


The maximum head coefficient obtained in these tests was about 3.6. 
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A chronological list of the test runs is given in Table IV. 
Also included are the important changes made on the TIR, the duration 
of each run, and any difficulties occuring during the run. The chief 
mechanical difficulties encountered were fluctuations in input pressure 
due to fluctuations of the axial compressor, slippage of the dynamometer 
torque arm during the run, and force capsule calibration changes due to 
changes in capsule temperature. These difficulties as well as several 
other minor difficulties were easily overcome by minor design changes. 
Calibration checks were made on the force capsules and the dynamometer 
arm poSition at the end of each run. With this information and the sta- 
bility of the input pressures during the runs, the author feels that 
runs 56 through 63, run 68 and run 78 yielded the most accurate perform- 


ance data. 


6. Results _and Discussion of Performance Tests 
General 

A complete set of raw data and reduced performance data for TTR 
test runs 51 through 80 is filed in the Turbo-Propulsion Laboratory 
Office, Reduced data for test runs 58 through 63, 68, 77, and 78 are 
contained in Appendix III. The performance parameters used in Appendix 
III have been defined in Sec. 4, In order to illustrate particular 
phenomena, some of the data from Appendix III havebeen graphically 
displayed in Figs. 18 through 56. Figures 18 through 29 are plotted 
to show the influence of axial clearance from the stator to the rotor, 
Ax, on the turbine performance. The influence on the turbine perform- 
ance of the radial clearance between the rotor and the shroud,Af, is 


illustrated in Figs. 30 through 36. Figures 37 through 46 are graphical 
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comparisons between the results of test runs 59 and 68, and perform- 
ance characteristics predicted for the Mod II Turbine by the three- 
dimensional prediction method described by Teeeiaoa 
Influence of axial élearance 

The most important finding of this study is that the total to 
static stage efficiency increases with increasing axial clearance up 
to an axial clearance,A x, of about 1.000 inch. This was found to be 
true at both radial clearances used in these tests. Figures 18 and 
19 show efficiency versus referred rotor speed at a stage pressure 


ratio, Pe [Pos of 1.3 for radial clearances, Af , of 0.033 inch and 


0,015 ee respectively. Data scatter in Fig. 18 precludes a quanti- 
tative analysis of the increase in efficiency with increase in axial 
clearance, But in Fig. 19, the data are smooth enough to state with 
certainty that the total to static stage efficiency increases about 
one point with an increase in axial clearance, Ax, from 0.200 to 
1,000 inch, Similar results are shown for a pressure ratio, (Po /P,)s 
Gt ernie ie see O and 21, . 
Equation (40) shows that for a given stage pressure ratio, Py IP 5, 
O 


a given inlet total temperature, T,. , and a given rotor speed, N, the 


L 
O 


stage efficiency is a function of the dynamometer torque, My and the 
mass flow rate, W, The influence of varied axial or radial clearance 
on efficiency can then be defined as a function of changes in mass flow 
rate and changes in dynamometer torque. By comparing referred mass 


flow rates, the effects of varied inlet temperature conditions need not 


be considered. Moreover, the changes of mass flow rate with changes in 





Heerison, R. G., An Analysis of Single Stage Axial-Flow. Turbine 
Performance Using Three-Dimensional Calculating Methods, (NPGS Thesis, 
Sept. 1967). 
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rotor speed are quite small. Therefore, for a given radial clearance, 
the influence of changes in axial clearance on the flow rate can be 
depicted by plotting referred mass flow rate as a function of pressure 
ratio. This was done for a radial clearance, Af, of 0.015 inch in 

Fig. 23 and for a radial clearance, AF, of 0.033 inch in Fig. 22. The 
conclusion drawn from Figs. 22 and 23 is that increased axial clearance 
has little effect on the mass flow rate. 

Similarly, to show the influence of varied axial clearance on the 
torque produced by the turbine, the referred dynamometer moment was 
plotted against the referred rotor speed. These plots are shown in 
Figs. 24 and 25, It can be concluded from these graphs that the increase 
in efficiency with an increase in axial clearance is primarily the result 
of an increase in produced torque. 

To summarize the fia three paragraphs, the Mod II Turbine will 
Operate most efficiently at an axial clearance, Ax, of 1.000 inch. 

At any other axial clearance the torque will decrease and the mass flow 
rate will remain about the same. 

During the test runs, the most noticeable effect as the axial 
clearance was increased was an increase in the static pressure at the 
stator tip and a decrease in the static pressure at the stator hub. 
This phenomenon is apparent in all test data. The pressure change 
shows up in the reduced performance parameters as a change in the 
theoretical degree of reaction. The theoretical degree of reaction 
increases at the tip and mean radii and decreases at the hub radius 
with increasing axial clearance. To show a typical result,the theo- 
retical degree of reaction at the tip and hub were plotted in Figs. 27 


and 28 for run 78. 


ge 


That the efficiency should increase with an increase in degree of 
reaction is quite consistent with the data shown by Fig. 26 which is 
taken from Teed The isentropic head coefficient for these tests 
varied from a kK. of 1.5 to about a ke of 3.5, and the theoretical 
degree of reaction at the mean was about 0.25. Referring to Fig. 26, 
an increase in a would be expected to increase the efficiency. The 
physical reason that the efficiency increases with an increase in degree 
of reaction is that the flow through the rotor becomes more accelerated. 
And a more accelerated flow means less opportunity for flow separation 
to occur. Knowing this, the phenomenon described above leads to the 
conclusion that losses are occurring in the Mod II rotor due to decelerated 
flow conditions. 

From 4x of 1.000 to 1.500 inch there was a further increase in a 
but a slight reduction in efficiency. This would indicate that, in this 
range of axial clearances, the increase in losses due to the increased 
length of the boundary layer along the shroud and/or other effects out- 
weighg the increase in efficiency due to the increase in the mean degree 
of reaction. 

The reason that the pressure distribution between the stator and 
rotor changes with a change in axial clearance is not entirely clear. 
One likely possibility is that the blunt leading edges of the rotor 
blades cause; an adjustment in the velocity field upstream of the rotor. 
Since the relative and absolute velocity fields between the blade rows 
are subsonic, the presence of the rotor can be felt upstream in the form 
of an induced pressure. If this is the phenomenon which raises the 


efficiency of the machine, the rotor must induce a pressure field which 


Ovavra, M., H., Aero-Thermodynamics and Flow in Turbo-Machines, 
New York, London: John Wiley and Sons, Inc., 1960, pg. 436. 


34 


guides the flow particles into the rotor flow channel at a more optimum 
angle as the axial clearance is increased. It has been demonstrated 

with stationary cascades on a water table that the axial distance between 
blade rows greatly influences the stream paths of the flow particles. 
However, the flow through a stationary row and a moving row of blades 
might produce different conditions. 

The change in pressure distribution may be due entirely to unsteady 
conditions caused by the moving rotor blades passing through the wakes 
of the stator blades. If this is the case, these phenomena will not be 
observed on a rectilinear cascade. In an unsteady flow,so-called 
Reynolds stresses are produced which must be added to the force field 
calculated from mean momentum flow. Experimental investigations about 
the magnitude of Reynolds stresses in turbomachines Seem to be non- 
existent in the available literature. Further studies with the TTR and 
Mod II Turbine should provide some answers to these important unsteady 
flow questions. 

Influence of radial clearance 

The total to static efficiency increased from one to four percent 
when the radial clearance between the rotor tips and the shroud was 
decreased from 0.033 inch to 0.015 inch. Figures 30 through 33 show 
the efficiency as a function of the isentropic head coefficient for 
pressure ratios of 1.3, 1.4, 1.5, and 1.6. There is no clear relation- 
ship between the efficiency change and pressure ratio, RPM, or axial 
clearance. It can be concluded, by comparing Figs. 18 through 21 and 
Figs. 30 through 33, that the efficiency of the Mod II Turbine will 
increase about three percent, if the radial rotor tip clearance is 
decreased from 0.033 inch to 0.015 inch. This increase is apparently 


independent of RPM, axial clearance and pressure ratio. 
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The most interesting finding was that the increase in efficiency 
with decrease in radial clearance is mostly due to a decrease in mass 
flow rate. That is, if all other conditions are the same, the Mod II 
Turbine will develop the same torque at a radial clearance of 0.015 
inch as it will at 0.033 inch but at less mass flow rate. This fact 
is graphically portrayed by Figs. 34, 35, and 36. Figure 34 shows the 
referred mass flow rate as a function of stage pressure ratio for both 
radial clearances of 0.015 and 0.033 inch. Figures 35 and 36 are plots 
of referred dynamometer torque as a function of referred RPM for various 
Stage pressure ratios at radial clearances of 0.015 and 0.033 inch, 
respectively. By comparing Fig. 35 with Fig. 36 it can be seen that the 
torque developed by the turbine is nearly independent of radial clearance. 
It should be noted that the pressure ratios portrayed in Fig. 35 are 
slightly different from those in Fig. 36. This fact must be considered 
on comparing the two figures since the torque developed by a turbine is 
very dependent on the stage pressure ratio. 

Figures 37 through 40 show the experimentally determined stage 
efficiencies on the same plots with predicted efficiencies as a function 
of referred rotor RPM. The predicted values are depicted by the curves 
and the experimental values by the plotted points. The experimentally 
determined efficiencies can be seen to be from one to three percent 
lower than the predicted values at both radial clearances, This dif- 
ference is mostly attributable to a difference in measured and predict- 
ed mass flow rates. Figures 31 and 42 are plots of the mass flow rates, 
and Figs. 43 and 44 are plots of the dynamometer torque. The predicted 
and measured torque data are generally quite closely in agreement. 
Likewise the predicted and measured power data, shown on Figs. 45 and 46, 


agree very well. 
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7. Description, Results and Discussion of Flow Surveys 


General 

Flow surveys were conducted with temperature and pressure probes 
before the stator, between the stator and rotor, and at the rotor 
discharge. The surveys were conducted to gain data from which compar- 
isons could be made with predicted flow properties and flow properties 
calculated from the mean streamline analysis method presented in 
Section 4. In addition, the flow survey ahead of the stator was 
conducted to gain a measure of the validity of the basic assumptions 


of the mean streamline analysis. 


Flow survey upstream of the stator 

A number of flow surveys were made with the two Kiel probes 
located upstream of the stator entrance. Total temperatures and total 
pressures were recorded at one-tenth inch increments between the radii 
of 3.7 inches and 5.0 inches from the TTR axis. Figure 47 depicts the 
results of a typical pressure survey with the Kiel probe on the left- 
hand side. The total pressure can be seen to vary about 2.5 percent 
between the stator hub and tip radii. The same distribution of pres- 
Sure was found by the survey data of the other Kiel probe. The tem- 
perature data from all traverses showed the total temperature was nearly 
constant at all radii. Therefore the assumption that the fluid prop- 
erties at the stator entrance are uniform is not entirely satisfied. 

The reason for the low total pressure at the hub radius is that 
the flow particles are accelerated around the hub curvature much more 
than the particles which follow a path at a greater radius. Since a 
constant loss coefficient is associated with the screens, the particle 


which passes the screen with the greatest velocity will incur’ the 
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greatest increase in entropy. Further, since no appreciable amount 
of energy is transferred in the process, all particles enter the 
stator at the same energy level (i.e.,total temperature). However 
the particles near the hub are at a higher entropy level and there- 
fore at a lower total pressure; thus, they have less energy available 
for the expansion across the stage. 

The non-uniform pressure causes two problems with the mean stream- 
line analysis. First, the total pressure measured ahead of the 
stator must be meaSured at such a location that it represents the 
pressure at the mass averaged mean streamline radius. Secondly, the 
loss coefficients determined by the analysis will not be exactly 
representative of the loss coefficients which would exist for uniform 
flow conditions. 

The positions of the fixed total pressure probes are at very nearly 
the same radius as the mass flow weighted mean streamline measured at 
the stator exit as shown in Fig. 47. It is therefore felt that the 
total pressure is very nearly representative of the mean condition. A 
design change of the TTR to allow for flow surveys immediately up- 
stream of the stator blades would permit the researcher to make a more 
quantitative analysis of the flow in this region. The one-dimensional 
analysis is designed to yield first-cut estimates of loss coefficients 
and other design parameters. Therefore, it is felt the non-uniform 
conditions at the stator entrance do not significantly detract from 
the usefulness of the calculated results. 

The assumption of axisymmetric flow conditions is apparently quite 
valid. The differences in pressure measured by the right and left Kiel 


probes at any radius were less than one percent of the absolute pressure. 
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Further, the maximum variation in pressure indicated by the six 
fixed probes was always about one percent. A typical result is 
shown in Fig. 48. 

Flow survey at the stator discharge 

Pressure and temperature surveys were conducted at the exit of 
the stator nozzles during runs 64, 65, 66, 77, and 80. Two 5-hole 
flow probes built by the United Sensor Corporation, type DA 125, were 
used on different occasions for the pressure surveys. These probes 
measure yaw and pitch angle, and total and static pressure. The total 
pressure and yaw angle are considered exact as read. However, the 
static pressure must be corrected for pitch angle, Mach number, and 
immersion effects, Additionally, the pitch angle must be corrected 
for Mach number and immersion effects. Calibration curves supplied by 
the manufacturer are filed in the Turbo-Propulsion Laboratory Office 
by instrument serial numbers. Since the Mach number of the flow at 
the stator discharge was much larger than the Mach number at which 
the probes were calibrated, it was necessary to linearly extrapolate 
the calibration data. The temperature surveys were conducted with a 
locally manufactured, shielded Iron-Constantan thermocouple. 

During run number 66,traverses were made at six peripheral positions 
from stator blade number four to blade number three. Ten data points 
were taken at each peripheral position. The position of each data point 
is shown in the plane of the stator blade trailing edges in Fig. 49. 
The calculated velocity distribution from the hub to the tip is shown 
in Fig. 50 for data points 21 through 30 and 31 through 40. These 
points were taken about half-way between the blades. The continuity 


equation was checked by 


he 
W = f, ery? Va, dr (52) 


oy, 


Equation (52) was integrated in an approximate manner by dividing 
the flow annulus into small increments, calculating the density and 
axial velocity that existed at the mean radius of each increment, and 


Summing as 


h i? 
, 7 4 
A G 14 
The only quantity on the right side of Eq. (52) which was neither 
known nor calculated from the measured flow probe data was the static 
temperature in the density relation. This temperature was found from 


the energy equation as 


2 
T, = h, 7 Wha Teg ae 

The results of this calculation are listed in Table V. The mass 
flow rate calculated from the survey data is 4.718 lb _/sec. The mass 
flow rate measured by the flow nozzle reduced by the labyrinth leakage 
flow is 4.430 lb _/sec. Thus the survey probe data apparently yield 
a mass flow rate which is about six percent too high. However, it 
must be kept in mind that the position of data points 21 through 30 
was in the region where the flow is nearly isentropic. Therefore it 
can be assumed that the mass averaged axial velocities are six percent 
smaller than the velocities listed in Table V. This reasoning can 
lead to the calculation of an approximate loss coefficient as follows. 


The mass averaged axial velocity is found as 


Rr 
a We _ (4.430) = 199 (fh /co, ) 
= Va Wy > 1 Ga.718) ladle (55) 
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View = axial velocity from Table VI at the radius, R which 


divides the flow channel into equal increments of mass flow. 


W = flow rate from Table VI. 





ox 
W = measured flow rate. 
Then 
Var 199 
yee eo = 585 (Fthec) (56) 
1 ~ Coe, Cos 70.1° 
where 
* 
Vi = mass averaged stator exit velocity 
x = measured stator exit angle at radius, Rg 


The stator loss coefficient can now be calculated as 





& 
“an Vi /eqJ op _ 4 BBS 723217) 18.6 _ 9. OG! 
s 


CC ie a _ 


Py = static pressure measured by the flow probe at RW 


Er = mean total pressure at the stator entrance 


| 
tl 


total temperature measured at the stator entrance. 


The temperature surveys at the stator exit indicated that the total 
temperature was very nearly constant at all radii. 

The results of the survey data are graphically compared to the 
predicted results on Figs. 55 and 56. The absolute exit angles, &,; 
coincide very closely with the exit angles computed in the mean stream- 
line analysis. However, the mean velocities from the mean streamline 


analysis were consistently higher than the measured velocities. 
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Flow survey at the rotor discharge 

Pressure and temperature surveys were conducted at the rotor 
discharge during run numbers 67 and 76. During run number 67 the 
total inlet temperature, qT, , total pressure, P. , and the RPM were 

O 
set to the same values as used in the traverse - the stator exit 
conducted during run 66. Therefore, the results from run 67 can be 
used with the results of run 66 to determine the actual velocity 
triangles at any blade radius. Similarly the inlet conditions of 
run 76 match run 75. The power developed by the machine can then be 
determined and compared to the power measured. Likewise, other 
performance variables can be computed and compared with the results 
of the mean streamline analysis and the predicted values. 

Figure 52 shows the velocity distribution at the rotor exit, 
measured during run 67, and computed by the method described previously 
for the DA-125 type probe. On the other hand, Fig. 53 shows the 
velocity distribution computed with the following assumptions: 

1. The total pressure, ew and total temperature, qT. , are 

assumed correct as measured by the probes. : 

2. The static pressure, Po» is taken as the atmospheric 


barometric pressure. 


3. The static temperature is then computed from 


ay 
bees (Typ * (58) 


4. The velocity is computed as 


V, = [2gSep (Tea -Ta>] 
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It is felt this method of data reduction is more accurate 
than the method which uses the probe calibration curves supplied by 
the vendor. Since immersion effects on static pressure measurements 
can be quite large , and are not taken into account,the calibration 
curves cannot be considered exact. 

Since the absolute rotor exit angle, ale has been measured, the 


axial and peripheral velocity components may be computed from 


V.= V2 Cos Xs ee, 9) 


and 


Ww. > Vv. Sin > (60) 


Also, the relative velocity components are obtained from 


Wa, = Varo (61) 
and 
Wy, = Vu, ~ V2 (62) 
Then, the relative discharge flow angle is 
B, = cintag lu (63) 
a2 


Table VI lists the results of run 67. The mass flow rate was 


computed as 
k 
W =a 2G Va, AP (64) 
he 


Because of the shape of the pressure probe, data points could be 
taken only at greater distances from the shroud than 0.16 inch. 

This condition made it necessary to assume the velocity distribution 
between the shroud and the meaSuring stations closest to it. Similar 


to the assumption used to develop Eq. (89), Appendix I, the axial 
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velocity distribution was considered to be linear from Zero at the 
shroud to the nearest measured velocity. With this assumption, the 
mass flow rate was determined to be 4.42 lb _/sec. This value coincides 
closely with the mass flow rate of 4.41 lb_/sec measured by the flow 
nozzle. 

The radius which divides the mass flow into two equal parts was 
found to be 4.361 inches. By comparison, the root mean square radius 
used in the computer program for the rotor is 4.187 inches. This 
difference causes the derived peripheral velocity from Eq. (24) to be 
in error by about six percent. Therefore, all quantities calculated 
by Eqs. (25) through (33) will be in error also. Because of this 
discrepancy, the rotor exit velocities, flow angles, and loss coef- 
ficients obtained by the mean streamline analysis method should be 
considered only as approximations of the actual values. 

The power produced by each increment of rotor area was computed 


from 


AHP = ato U (Vu, - Woz) G65) 


The formulas for calculating the variables in Eq. (65) are listed in 
Table VI. 

From a summation of these increments, the total power developed 
by the rotor was determined to be 63.3 horsepower. The power calcu- 
lated by Eq. (38) from the measured dynamometer torque and the RPM was 
62.1 HP for the data from run 66 and 63.1 HP for the data from run 67. 

The results of the flow survey from run 67 are compared to the 
predicted results given by Haneteen on Figs. 55 and 56. It is felt 


that the discrepancy between predicted and measured exit angles is 


7 : ; 
Haeersen, 1Oc., Cit, 
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caused by the separated flow regions in the rotor which were not taken 
into account in the prediction method. In Fig. 5 dark areas can be 

seen at the tip and hub of the rotor blades. These areas consist of 

oil and dust which accumulated on the blades during the runs because of 
separations at these locations. Also in Fig. 57 dark areas can be seen 
along the blade near the hub, as well as dark lines that point radially 
outward. Apparently thick boundary layers or separated flow regions 
exist in these areas. The dark lines are probably caused by dust 


particles being thrown outward by centrifugal forces. 


8. Recommendations 

The following recommendations are not presented as original 
thoughts of the author, but rather as a listing of the logical steps 
which should be undertaken as a continuation of this study. 

1. Blade shapes geometrically similar to the Mod II blading 
should be tested in the Rectilinear Cascade Test Rig of the Naval 
Postgraduate School. The loss coefficients found in the present 
study could be validated and the interesting changes of the pressure 
distribution with axial clearance could be investigated. 

2. The non-steady flow conditions caused by the rotor passing 
through blade wakes of the stator should be examined in greater detail. 
The wake patterns could be measured by hot wire anemometers now avail- 
able. The change of the pressure distribution with axial clearance 
may be related to non-steady flow phenomena. 

3. The TTR should be modified to obtain more uniform flow con- 
ditions at the stator entrance. A fix that should be attempted is to 


remove the turbulence reduction screens and to insert anouter wall 
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contour in the channel leading to the stator entrance, If this change 
is unsuccessful in reducing turbulence and increasing uniformity, then 
the TTR stator assembly should be fitted with a larger diameter inlet 
pipe. Also, it would be very useful to have the capability of conduct- 
ing a flow survey immediately upstream of the stator entrance. 

4. Further performance tests should be conducted on the Mod II 
Turbine. Particularly, tests at higher head coefficients should be 
performed after the water brake dynamometer becomes available. It may 
be determined that the relationship between efficiency and axial clear- 


ance is a function of head coefficient. 
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TABLE I 


MOD II TURBINE DESIGN PARAMETERS 





LTEM SYMBOL UNITS STATOR ROTOR 
Number of Blades 
Blade Height | h 
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TABLE II 


MEASURED DATA FOR TURBINE PERFORMANCE TESTS 


TRANSONIC TURBINE TEST RIG 


Desired Measured 
OQuantit Item 
FLOW RATE 


Flow nozzle total pressure 


Flow nozzle total temperaturé 


Pressure differential across 
flow nozzle 


Turbine plenum total pres- 
sure 


STATOR ENTRANCE PROPERTIES 


Stator entrance total pres- 
sure 


Stator entrance total tem- 
perature 


STATOR DISCHARGE PROPERTIES 
Stator tip static pressure 
Stator hub static pressure 


Force exerted on stator 
assembly by the axial force 
capsule 


Moment exerted on stator 
assembly by the stator 
torque capsule 


Static pressures at the end 
of the shroud 


Axial force exerted on the 
stator assembly by the 
closure plate 


Moment exerted on the stator 


‘Measured 
Units 


in. Hg. 
millivolts 
in. H,0 


in. Hg. 


in. Hey: 


millivolts 


in. Hg. 
ine Hg ° 


counts 


counts 


| in. Hg. 


counts 


counts 


assembly by the closure plate 


ROTOR DISCHARGE PROPERTIES 


Static pressure in the test 
hood 


Dynamometer torque 


Rotor rotational speed 
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ine He. 


counts 


RPM 


Symbol 


Computer 
Character 


PNDZ 
TNDZ 
DH 


PSPL 


PPPL 


TTPL 


BLIP 
PHUB 
FAX 


TORQ 


yl hora) 74 8 


CLFAX 


CLTRQ 


PHD 


DYNAR 
RPM 


TABLE III 


MOD II TURBINE TEST RUNS 1967 





JAN.-AUG. 
E 
Run No. ar Ax to/, N 3 Test 
(in) (in) 2 (RPM x 10 ) Purpose 
(approximate) 
o 51 0.033 0.410 ee 13-20 + 
BE 32 0.620 1.4,1.5 12) - 
* 53 0.620 oe ae 11-18 ati 
* 54 O74 LO Oe 15-18 + 
ee 2) 3 0.410 0 I 
* 56 0.410 Li 13-19 il 
* 57 0.4150 1.4 13-19 2 
* 58 0.620 A i ls 10-17 
* 59 1.000 1.3,1.4 ime. 6 LO-17 
* 60 0-200 1.3, le@iler, Peo 10-17 


208 0.410 
* 62 . 0.410, 
* 63 Le YC, 


mle te 





3 + 

3 + 

3 + 

3 + 

5 + 

3 + 
o 64>) 0201 1.000 1.4 14.6 3 
o 65 1.4 14.6 5 
o 66 1.4 14.6 3 
o 67 1.4 14.6 4 
* 68 eBid idesub, Dd.) 6 11-18 + 
* 69 0.410 re ee 11-19 + 
m7 0.410 is 14-19 + 
= 0.200 es, Le 11-17 + 
* 72 0.410 L . Saeed 11-18 + 
aes: O2Z200 or, ee sf 
o 74 02208 deus chat = 
Ove 0.410 Lea 11-5 5 
o 76 0.410 pigs a Pe 11-19 4 
o 7/7 0.410 & 1.000 L Sales 11-19 + 
o 78 2008. Gies e00 er. 3 11-18 2 
wee 200 he oul ee 11-17 + 
o 80 410 lbs, 10-17 a 
Oo = without exhauster P= poy 1 = Dynamometer check 
* = with exhauster 2 = Spring capsule dyna. 
+ = performance characteristics ee x 3 = Stator exit survey 
- = hood check 4 = Rotor exit survey 


nO 


TABLE IV 


CHRONOLOGICAL RECORD OF MOD II TURBINE TESTS 1967 


Run No. Date Duration Remarks 
(Hrs ) 

or 1-26 SIMS) namometer arm slipped, put set screw in arm 
32 2-3 25D Tested hood for leaks, tested dynamometer 

2-9 Tested waterbrake, slight bearing problem 
By: 4-3 oie) Dyna. capsule 25 counts off at shut down 
54 4-24 Ss. erp = 1.6 & 1.8, dyna. capsule off 30 counts 
55 4-25 Checked dyna. capsule at various temperatures 


by heating with a lamp. - very temp. sensitive 
-50 to +170 counts 

Installed constant temp. environment system 
around dynamometer capsule 


56 5-11 S510 Performance test at Resp, = 1.4 with dyna. 
capsule; took tare readings during run. Dyna. 


checked good at shut-down 














57 5-16 25 Used spring capsule to check dyna. readings 
ob tained during run 56 - satisfactory 
28 3t 6 Good data - no apparent problems 
59 6-1 6.0 A " " T T 
60 6-5 oD au " " " " 
61 6-7 a0 a " " " " 
62 6-9 4.0 i UL i" i" " 
63 6-13 5.0 " " i" " '" 

toe Replaced shroud with new shrou@$ Ar = 0.015 in.; 


cleaned all probes & lines with high press. N,. 
Unable to calibrate closure plate axial force 
device. Designed new device & assembled mach 


w/o closure plate instrumentation 

64 6-23 6.5 Temp. & press. traverse at stator exit. Y.C. 
DA125_ S/N 926 

65 6-26 3.5 11 1 1" 1" 11 1! 1 





66 6-27 Bree, 60 data point pressure traverse at stator exit 
with DA125 S/N 926, 10 data points with S/N 928 


for comparison. 
67 6-28 a0 Temp. & press. traverse at rotor discharge. 
Y.C. D.A. 125, S/N 926. Two peripheral positions 
checked about 30° apart.-readings less than 
+ 0.25 cm. HO. difference 
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TABLE IV (cont.) 


68 7-5 oe) Good performance data - no apparent problems 
69-76 7-10 44.0 All dyna. readings during these runs are 
to possibly in error. Set screw had worn a grove 
8-5 in the shaft causing slippage. Discovered 
during run #77. Also A.C. compressor very 
unstable. 
7-20 Installed new closure plate axial force 
flexture, Calibrates good, 
AS, 7-31 Traverse before stator with left & right Kiel 
probe, Transverse aft of stator with S/N 928, 
76 Traverse aft of rotor with S/N 928 
19 Discovered torque arm movement during run - 


made calibration of movements. 


78 opring capsule dynamometer used 





79 Checked capsule calibrations 
80 Traverse aft of stator with S/N 928 
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FIGURE 1 


MOD II TURBINE STATOR ( After tests) 


(View showing stator entrance) 
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FIGURE 5 


MOD II TURBINE RCTOR (After tests) 


(View shov ing rotor exit) 
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FIGURE / 


MOD II TURBINE BLADE PROFILES 
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FIGURE 11 
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FIGURE 15 THERMODYNAMIC PROCESS OF FLUID IN AN 
_____ AXIAL TURBINE STAGE | 


STATOR EXIT PLANE 


> as Vu, ROTOR EXIT PLANE 








FIGURE 16 VELOCITY DIAGRAM OF A TURBINE STAGE 





FIGURE (7 


VECTOR SYSTEM FOR THE ROTOR ASSEMBLY 
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FIGURE 45 


99 





3 I¢ 15 6 TT 18 


x 1073 (RPM) 
REFERRED POWER VS REFERRED RPM 
Ar= 0.018 In. Ax 21,000 In. 
FIGURE 46 





$1X3 204035 
iy 9U1}W0813S 
udeyy PanOiem 
0) 4 S8OW 


S 9 


ei= 94 gzeNNY 3g0Ud T3IN 1437 HLIM 
JONVYLNZ YOLVLS 1V ASAYNS 3YNSS3¥d 


Ze 3YNDSDI4 


a’ 








Po) +. 


{ 
9q04g 


“p 
9INSS81Ig [040] pexi4 J ov 





| 
nf @snsseig uDEW lob 


TP peiubiem 013 sso 


GE =" 


J! 
: / 
9S . ] 


/ St 


+ 


] 
S 
OPS NE SENN 
(6yyu! ) ‘, 8INSSO1g |DJO] 


Radius (inches) 


101 


8b 3yYNSIS 


Giz NNY S380Ud Q3Xi4 3HL AG 
Q3UNSV3SW SV JONVYINS YOIVIS 3HL LV 
NOILISOd TVYESHdIdad HLIM NOILVIYWA 3YNSS3Yd WWLOL 


UOINISOd |048YdI28qy 403DkS SuUIquN) || POW 


o09E e00E o0ve 008! 00d! 009 me) 


do) do} 
£1 o> te 8¢ 


~~ = = = 
oo o— 
=_ = 


-_=a=a= ae 
= * 
= 


vA 
P 
(6y'u))°y esnssesg jojoL 


O 
Sg 


woe8sjsdp 
Buiyoo7 MBIA 





102 


PRESSURE PROBE POSITIONING AT MEAN RADIUS 
RUN “66 
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L.E. 
Blade*3 


STATOR DISCHARGE PRESSURE SURVEY RUN#66 


(Flow Probe DA—|l25 *926) 
FIGURE 49 
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FIGURE 5l 
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ABSOLUTE FLOW OUTLET ANGLES 
AS FUNCTION OF RADIUS 


(Ax=1.000in., Ar=O0.015in, P /P=140, RPM/ 16 =13,934) 


FIGURE 55 
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APPENDIX I 


FORMULA DEVELOPMENT FOR MEAN STREAMLINE ANALYSIS 

A. The peripheral component of the velocity at the stator exit is 
determined by applying the moment of momentum equation to the fluid 
contained in the stator assembly. The assumed control volume sur- 
rounds the fluid in the assembly. It is assumed that significant 
flow separations do not exist, and the fluid shear stresses at the 
entrance and exit of the control volume are thus negligible (see 
jes) With this assumption the moment of momentum equation may 


be written for this application as 


M=-f F xYdm,, + f,, Bx Mame, (66) 


where 


M = Moment applied to the control volume by the stator assembly. 
- = Surface where fluid enters the control volume. 
Ss) = Surface where fluid leaves the control volume. 
The vector conventions for Eq. (66) are depicted on Fig. 13. Since 
net pressures which could produce moments on the stator assembly 


do not exist, pressure integrals have been omitted from Eq. (66). 


In the axial direction the component of M is 


M; = -F f ae Vis dm,, + fe r Vu, ding, (67) 





Sumare: M. H., Aero-Thermodynamics_ and Flow in Turbomachines, New York, 
London: John Wiley and Sons. Inc., 1960, p94. 
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The air enters the stator assembly from the plenum in a radial 
direction. Therefore the peripheral velocity component, Vy , is 


oO 
zero. Thus 


Since the average flow conditions are assumed to exist at the 


mean radius, Eq. (67) can be integrated to yield 


M: = Rm Vu, m he, 
The moment applied to the inside of the shroud by the rotating 
fluid downstream of the stator exit is considered small. Thus, the 
moment applied to the fluid by the stator is the sum of the moment 
applied to the stator assembly by the stator assembly moment capsule, 
and the moment applied to the stator assembly by the closure plate} 
hence 
Mi = Mer Me. (69) 
where 
M, = Moment applied to stator assembly by the stator assembly 
moment capsule. (ft-1b,) 
M.. = Moment applied to stator assembly by the closure plate. 


CL 
Substituting Eq. (69) into Eq. (68) and with m =W/g, there is 


Ve = (M+ Med (70) 

UI = WW Ren 

B. The axial velocity component at the stator exit is obtained by 
applying the conservation of momentum law to the same control volume 


as used above. For this application the momentum equation may be 


written as 


(a = Ve dm,. + Ie V,drh., + SARS. - SA RAS, Gi) 


Ba 


where 

F = Force applied to the fluid by the stator assembly. The 
vector field conventions are shown in: Fig. 13. 

Since neither fT, nor V, have a component in axial direction, the 


vector component equation in the axial direction is 


= tf, Va, dirs, - f,, RAS (72) 


The pressure integral cannot be integrated unless the pressure 
distribution at the stator exit is known. The pressures at the tip 
and hub of the stator are recorded and the distribution is assumed 
to be linear. The static pressure at any radius, r, at the stator 


exit is therefore assumed to be 


=f, + (Li *)(R- BR) (73) 


where 
P, * Static pressure at radius, r, at the stator exit. 
= Static pressure at the stator hub, 
a Static pressure at the stator tip. 
= Radius at the stator hub. 
r, = Radius at the stator tip. 
The pressure integral in Eq. (72) can now be integrated as 
fRds, sam[ Rip (hn) + Rm (Rot (74) 


= KRW Ge )) 
oh 
and the velocity integral as 


if Va,dm, = Va, M = Va, Wa (75) 
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The axial force applied to the fluid in the stator assembly is 
reacted to by the sum of the forces applied to the assembly by the 
axial force capsule, the closure plate, and the net pressures; 
that is, 

Rak +k + 2, 22 
where 


F = Force applied to the stator assembly by the axial force 


capsule (1b ,) 


For = Force applied to the stator assembly by the closure 
plate (1b ¢) 
ae = Force applied to the stator assembly by the net pressure 


force (1b ,) 


Figure 13 shows the applicable areas used in Eq. (77) for 


determining the net pressure force. 


2A = ~ Prog An* te Ab+RAc+R Aet tr An _— 
With Eqs. (74), (75), and (76); Eq. (72) is solved as 
Vo, =F +F + DA - 27 ty, (ht ho )+ R-R 
re = Fn 


(78) 
nr? 73 z 3 Sho. 
Cee = Be - Ie yy WW 


The thermodynamic properties of the fluid at the stator discharge 
and the performance parameters for the stator can now be calculated. 


Referring to Fig. 15, the static temperature is 


2 
= Gi 5 v7 = Va, + Wy, (79) 
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From the equation or state the density is 


eC = T/ RT, (80) 


The equation of continuity can now be checked as 


W=CAM,, 5 A, =Aax Keg, (81) 


where 
Ay = Effective axial flow area at the stator exit. 
Ax = Actual axial flow area at the stator exit. 
l 
Kee = Factor to account for finite thickness of the stator 
1 
blade trailing edges. 
ae is given by Ten from cascade data as 
l 
4-27 (168 ty, ee ays (82) 
ee 102 ' ' 


where the variables in Eq. (82) are shown in Fig. 7. 

The axial velocity component at the stator exit may be obtained 
independently from the momentum methods mentioned previously by 
combining Eqs. (79), (80), and (81). First,the static temperature 
is obtained by combining Eqs. (79), (80), and (81). The resulting 


quadratic equation gives 


_ Joo KR 2WwiRe 7 Wi 1 ah 
= R? W* (UI -SAae Care “WI 1} (ee) 


Using Eq. (79) the axial velocity component is found as 


| 
Va, = reglep(h,-h)-W, J] * (84) 


eaten M. H., Unpublished notes for course Ae 432, Naval Postgraduate 
School, 1966. 
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C, The peripheral component of the relative velocity at the rotor 
discharge is found by the application of the moment of momentum 

equation to the fluid contained between the rotor blades. With the 
same assumptions as used previously, the component of the relative 


moment of momentum equation in the axial direction is 


M oe. Tile CR, Wu, + @Rt)drng + Ff, (Ra Wu, WRi Ome, (85) 


_ = 
where 

M, = Moment exerted on the fluid between the rotor entrance 

and exit. 

The vector field conventions are shown in Fig. 17. Once again, the 
shear stresses at the entrance and exit surfaces are considered 
negligible. Therefore Eq. (85) is only valid when large separated 
flow conditions do not exist. 


Substituting W/g for m and integrating 


M. = (Wu, 7 W.,) Rm Wo ay 
M, is a negative quantity and is equal to the moment applied to 
the rotor shaft by the dynamometer, My? plus the moment absorbed by 
the rotor and dynamometer bearings. The bearing losses were calculated 
and found to be of the order of one-tenth of one percent of the power 
produced by the Mod II Turbine. Therefore, the bearing losses were 


ignored. Substituting M) for (-)M, in Eq. (86) gives 


Wu. = Wu, - 9Mo/p wy (87) 


The relative axial component, W is determined by satisfying 


a2’ 


the equations of continuity and energy. The continuity equation is 


W= A, Wa, 3 Art Aax Kten (88) 
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A, = Effective axial area at the rotor exit. 


> 
It 


Actual axial area at the rotor exit. 


aw 
ul 


Factor to account for finite thickness of the rotor 


blades at the trailing edges. 


Equation (82) must be modified slightly for the rotor because 
of rotor tip clearance effects. This is accomplished by assuming one- 
half of the area between the blade tips and the shroud to be the effect- 


ive flow area of the tip clearance flow. ae is then 
2 


Ay 


3.5 @ Arh 42 (89) 
K,, 24 1- 23 (100% Y, 6 St a ooo 
1 3 = 
e. ras C 2) hh)? 
where t , S,, and a, are shown in Fig. 7 and where: 
ey 2 2 


Ar = Radial clearance between the rotor tips and the shroud. 


K 
Il 


Tip radius of rotor blades. 


—=sHuberadius of rotor blades. 


a) 
{ 


For the relative flow field, with no change of mean radius across 


the rotor, the energy equation can be written as 


2 z 
ir = iP te iets = iG + W359 ep (90) 
where 
Ty = Equivalent temperature (see Fig. 15). 
T, = Static temperature at the rotor exit. 
Wo = Relative velocity at the rotor exit. 


Rearranging Eq. (90) and replacing W, by the components W and W 


2 a2 u2> 


T,shewe  -wi,. - wi (91) 
gle "Sg Jep 7297 ep 
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where the unknowns are T. and W .. 
2 a2 


By the equation of state 


(92) 


GO AN 
where 
Po = Static pressure at the rotor exit = pressure in the test hood. 


Combining Eqs. (88), (91), and (92) and solving for T, by the quadratic 


Z 
equation 
_ gJep Ar $2 f 2w?Re Woe _ 93 
ee [1- BT ep Re Wad ep q- rere J} } oe 


Then by Eqs. (88) and (92) 


Wa. = W.Re Te (94) 


gle. 


APPENDIX II 


FORTRAN IV COMPUTER PROGRAM FOR PERFORMANCE DATA REDUCTION 
A. Inputs 

The input variables and their units are explained in comment 
statements in the executive program. The FORTRAN input statements are 
contained in Subroutine INPUT. 
B. Outputs 

Three output forms are included in the program. 


1. A printout of the input variables is accomplished by Subroutine 


2. A printout of the general performance results is accomplished 
by Subroutine OUTPUTA. 

3. A printout of the referred stage performance parameters is 
accomplished by Subroutine OUTPUT... 


Samples of the last two output forms are contained in Appendix III. 
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APPENDIX III 


REDUCED PERFORMANCE DATA 

A, Order of listings 

The referred performance parameters and general flow results for 
runs 58 through 63, 68, 77, and 78 are listed in numerical order. The 
data point number is Listed on the far left hand side of each even 
numbered page. Read horizontally across both the odd and even numbered 
pages for the results pertaining to that data point. 
B. Definition of parameters 

The listed parameters are defined in Sec. 4. 
Cc, Units 


Ll. The units of referred parameters are listed at the top of each 


2. All velocities are feet per second. 

3. All temperatures are degrees Rankine. 

4, All angles are degrees. Positive angles are measured in the 
direction of the rotor rotation vector. 


5. All other parameters are dimensionless. 
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TURBINE TYPE MOD II CONFIGURATION RADIAL ROTOR 1 
TEST RUN-NO. 77 DATE OF TEST 
POINT PRESSURE  ISENTROPIC EFFICIENCY REFERRED 
RATIO HEAD COREE BOI-STATIC™ Clow RATE 
(R=4.125 IN.) PERCENT LBM/SEC 
l 1.2988 3.0431 80.88 3.1039 
2 1.3011 2.5106 79.16 3.0898 
3 1.3008 1.8738 75.68 3.0719 
4 1.2991 1.3709 63.49 3.0765 
FOR POINTS 1 TU 4 AVG. PRESSURE RATIO= 1.2999 » MAX.D 
5 12400 2.9253 80.13 3.3456 
2 1.3986 2.5898 78.78 3.3330 
7 1.3992 2.1042 Ts 2T 3.3292 
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14 1.4012 - 5888 81.02 3.3410 
15 1.3982 2.1028 79.40 3.3254 
16 1.3989 1.7436 73.69 3.3250 
FOR POINTS 13 TO 16 AVG. PRESSURE RATIO= 1.3998 » MAX.C 
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